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As a prototype of systems bearing a localization-delocalization transition, the strand-separation 
(melting) process in a double-stranded biopolymer is studied by a mapping to a quantum-mechanical 
problem with short-ranged potentials. Both the bounded and the extensive eigenmodes of the 
corresponding Schrodinger equation are considered and exact expressions for the configurational 
partition function and free energy are obtained. The force-induced melting is a first order phase 
transition process, while the thermal melting is found to be second order. Some scaling exponents 
governing thermal melting are given. 
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DNA melting, the strand-separation of a DNA double- 
helix, is an issue with both practical biological signifi- 
cance (because it is closely related to DNA replication 
and gene transcription) and pure academic interest. The 
first realistic model of DNA melting was proposed by 
de Gennes in 1969 ^ and reintroduced by Peyrard and 
Bishop twenty years later Q. The essential advance of 
the de Gennes-Peyrard-Bishop approach, distinguishing 
it from earlier efforts based on Ising-like models is 
that (i) the continuous degrees of freedom for the con- 
figurational fluctuation of a DNA and (ii) the short- 
ranged hydrogen bonding between its two complemen- 
tary strands have been explicitly incorporated. The DNA 
molecule is considered as consisting of two flexible Gaus- 
sian chains with short-ranged on-site interactions be- 
tween them; its statistical property can be studied by a 
mapping to a weakly-bounded quantum mechanical prob- 
lem 1^,^. The corresponding Schrodinger equation has 
a flnite number of bounded eigenstates and a continu- 
ous series of extensive ones, representing respectively lo- 
calized and delocalized eigenmodes of structural fluctua- 
tions. In fact, it has been well recognized that problems 
such DNA melting, flux line depinning in superconduct- 
ing shells B, adsorption of polymeric materials onto a 
surface [0J|], and some wetting phenomena 0, are all 
governed by the competition between enthalpy-favoring 
localized states and entropy-favoring delocalized ones. 

We study the melting of a double-stranded biopolymer 
as a prototype of such many phenomena. Earlier stud- 
ies have focused mainly on thermal melting ||l|,^,p|-p^ , 
here both the thermal and the force-induced melting pro- 
cesses are discussed; and different some previous efforts 
P,^ HTl|JT^ which considered only the localized ground- 
state, we explicitly incorporate all the bounded and de- 
localized eigenmodes of fluctuations (we are grateful to 
Prof. D. R. Nelson for informing us that in ||l3| a simi- 
lar treatment has been performed in studying votex pin- 
ning). This makes it possible for us to obtain exact ex- 
pressions for the partition function and free energy of 
the system. We rigorously show that the stretch-induced 
melting is a first order structural phase transition while 



that induced by heating is of second order. The phase 
diagram for the double-stranded polymer is exactly ob- 
tained, and some scaling laws governing thermal melt- 
ing are given. The effects of sequence heterogeneity and 
base-pair stacking are also briefly discussed. 

The model double-stranded biopolymer is formed by 
two Gaussian chains. Each strand contains TV -I- 1 beads, 
with a harmonic attraction between any two consecutive 
ones; and between each pair of beads of the two strands 
with the same index there is a short-ranged interaction 
potential V . The model energy is 
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= J2 [<1^n - Rn-if + J(r„ - r„_i)2] + ^ y(r„). (1) 



Here, rl*-* denotes the position vector of the nth bead 



in the ith strand; and r,. 



and R,i 



(ri^-* -I- rn')/2 are, respectively, the distance between a 
pair of beads of the two strands and the ccnter-of-mass 
position of these two beads. For computational simplic- 
ity, in the following we discuss only the one-dimensional 
case of model Eq. (|l|). The principal conclusions of this 
work is independent of dimensionality as well as the par- 
ticular forms for the short-ranged potential V, since the 
underlining physics, the competition between enthalpy in 
the localized states and entropy in the extended states, 
is reserved. We flrst discuss the situation of symmet- 
ric potentials and assume the short-ranged attraction to 
be (5- form, V{r) = —^6{r). At index n = the two 
strands intersect each other, i.e., Rq — rg — 0. The 
statistical weight for the center-of-mass position at the 
other end to be equal to R is easily known: Zfj {R, N) = 
{(3K/'!TN)^/^e^p{-l3KR^/N), where f3 = l/fcsT with T 
being the temperature. 

The statistical weight Zr{r, N) for the relative distance 
r between the two strands at index N is governed by the 
Schrodinger equation ||l|,0] 



.(2) 



1 



dZr{r, N) 
dN 



PV{r) 



Zrir,N), 



(2) 



with the initial condition Zr{r,0) — 5{r). Equation ^ 
corresponds to a weakly-bounded statistical system with 
only a finite number of localized eigenstates (in the case 
of (5- form attraction used here, this number is unity); the 
ground eigenfunction of Eq. (||) is (j)b{r) oc exp(— 77|r|/2), 
where 77 = The statistical weight Zr{r,N) is ex- 

pressed as 



+ J dXexp[~-Ne{X)](t){X,r)(l3*{\,0), (3) 

— 00 

where 6(A) and 0(A, r) are the eigenvalue and the 




rj Nrf' . 



TT , I3k 2 

77exp( r 



and the total partition function is thus 

S(/,iV) =exp(7V(/3/)V4/3«) 



N' 
4^ 



eigenfunction of the extensive eigenstate of Eq. (g) 
with wave number A [^. For the force-fixed 
ensemble with an external stretching / acting on 
the A^th bead of the first strand, the partition 
function is E{f,N) ~ J dRjy J dr^ exp{(3fRN + 
{l/2)PfrN)ZRiRN,N)Zr{rN,N). Since the first term 
of Zr{r, N) is proportional to exp(— 77|r|/2), it seems that 
when f > 7]/ (3 this integral will turn to be divergent. Ear- 
lier studies considered only the ground eigenstate 
and therefore took such a divergence as signifying the 
occurrence of force-induced melting process. Actually, 
however, there is no divergence problem. After taking 
into account of the second term in Eq. (||) , this term is 
canceled out by a term resulted from the integral. The 
correct form of the statistical weight is 
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where ^0 = (2/v^) dt^ exp(-i2) j^^^^^ dy exp[-{y^- 
0^ Pt^ /if)] is an small quantity. Equation (H) shows that 
in the thermodynamic limit, the free energy linear density 
is g{f) = lim -kBTE{f,N)/N = -K-f^l3^/A - P/Ak 

for f < fc and -P/2k for / > /c, with = ^7/?. 
Therefore a first order phase transition occurs at the 
threshold force fc- The inter-beads distance for the first 
strand is f /2k for f < fc and //k for f > fc and a dis- 
continuity appears at fc. Similarly, the average distance 
between the end beads of the two strands is approxi- 
mately zero for f < fc and proportional to N for f > fc 
[rN = Nf/n, see Fig. 1). At fc, f^v = S/r; + N7j/2kP. 

The extension-fixed ensemble may be more directly re- 
lated with actual experiments it is much easier 
for one to fix the total extension of the first strand than 
to fix the external force. The statistical weight for the 
end-to-end distance of the first strand to be fixed at aN 
is equal to J dR J drd{aN - R - r/2)ZR{R, N)Zr{r, N). 
a is the inter-beads distance of the first strand. Based 
on Eq. (^), we find the free energy density to be g{a) = 
KCT^ — rf /Ak.(3'^ for a < ac, g{cr) = 'qa/ (3 — rf I2k.I3 for 
fc < cr < 2ctc, and g{<j) — ko^ 12 for a > 2ac, with 
"■c = 7/3/2. The free energy function is hence a piecewise 
smooth function. Correspondingly, the average force is 
/((t) — 2k(t for a < CTc, f(cr) ~ k'jP for CTc < cr < 2ac, 
and /(ct) = na for a > 2ctc (see Fig. 1). The occur- 
rence of a force platform may be striking . It indi- 



cates that as the inter-beads distance in the first strand 
reaches CTc, melting of the double-stranded polymer origi- 
nates from the end point (index N) and progresses along 
the chain until the whole polymer becomes separated. 
Stretch-induced melting can be termed as directional 
melting 18 1. The phase diagram of this system is 
shown in Fig. 1, it includes a double-stranded native re- 
gion, a single-stranded denatured region, and a coexist- 
ing region. This transition is caused by enthalpy-entropy 
competition, different from that discussed in Ref. [ p^ , 
which is caused by the appearance of two minima in the 
ground-state eigenfunction. 

The above model with symmetric potential does not 
exhibit thermal melting behavior. In the following we 
improve our model by changing the attractive potential 
in Eq. (|l|) into the following asymmetric form V{x) = 00 
for X < and V{x) = — jS{x — a) for a; > 0, where a is a 
characteristic distance. For such an asymmetric system, 
when the temperature becomes high enough, the local- 
ized eigenstate disappears. Hence it might be possible 
to qualitatively describe the thermal melting of double- 
stranded biopolymers. 

We focus on how the distance between the two strands 
changes with external stretching or temperature. It is 
convenient for us to assume that ro = a. For this revised 
model system, we find that the statistical weight for the 
relative distance is 
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In the above expression, r = aKjP'^ and C is the largest 
solution of C = 1— exp(— Ct). This equation has a nonzero 
solution only when t > Tc — I. When r < Tc the solution 
is C = 0. Since r decreases as temperature increases, the 
polymer's structure might undergo a transition at tem- 
perature T — T„i — yfRFfa/kB- 

When an external force is acting on the first 
strand, the total partition function is S(/, N) = 
exp(Af/3VV4K/3) /o°°rfrexp(/3/r/2)Z^(r,iV). This inte- 
gral is always convergent for any value of /, although 
the first term of Eq. (|^) scales as exp(— CTr/2a). The 
integrand in the second term of Eq. (||) has two poles at 
A = ±<^Tz/2, hence it will generate a term which precisely 
cancels out the first term in this equation. To obtain the 
analytical expression for the partition function N) 
is nevertheless a demanding task. We evaluate alterna- 
tively its Laplace transform: 

oo 

£[S(/,iV)](s) = / diVexp(-siV)S(/,Ar) 



8v/s(K/3)3a[exp(/3/a/2) - exp(-i 
[4sK^ - /32/^] [2^/^a - r(l - exp(- 



(7) 



The largest solution of the equation Af)](s) = 

corresponds to the linear free energy density of the poly- 
mer system. When the temperature is less than T^, the 
free energy density g{f) = —{(T)^/4:K(3'^a'^ — /^/4k for 
f <fc and g{f) = -P/2k for / > /c, where fc = Cr/ap. 
Thus, the external force will induce a first order phase 
transition at the threshold force fc, which decreases as 
the temperature increases (see Fig. 2). This is similar 
with what we have attained with the earlier model. At 
T = Tm, the threshold force decreases to zero. The sta- 
tistical behavior of the extension-fixed ensemble is also 
similar with that of the earlier model and the phase dia- 
gram is shown in Fig. 2. 

When there is no external force, the free energy density 
is -{Crf /AnP'^a'^ for T < and zero for T > r,„. 
At Tm , the free energy and its first order derivative with 
temperature is continuous but the second order deriva- 
tive is not, indicating that the thermal melting at T„i is a 
second order continuous phase transition, with a discon- 
tinuity in the specific heat. The order parameter for the 
thermal melting can chose to be the probability Pioc(n) 
for the distance of a pair of beads (with index n) of the 
double-stranded polymer to be less than the character- 
istic length a. For the thermal melting process, we can 
predict based on Eq. (0) that, as the melting tempera- 



ture T,n is approached from below, P/oc ~ {Tm — T)^ , 
with the critical exponent /3 = 1 |^^; it is also easy to 
obtain that as the temperature approaches T^, the cor- 
relation "length" in Pioc{n) between different indices n 
and 7i' scales as \Tm — T\~'^, with the critical exponent 
7 = 2 !§. At T„, {Pioc{n)Pioc{n'))c - - 
|pO| , but the critical exponent fj is difficult to be obtained 
by the present asymmetric model. Nevertheless, we no- 
tice that the phase diagram for the symmetric (Fig. 1) 
and the asymmetric model (Fig. 2) is identical albeit 
that the symmetric model has = oo. Therefore, it 
should be possible for us to obtain a good estimation of 
77 based on the symmetric model by artificially assuming 
7 = jo{Tjn — T)/Tm for T < Tm and 7 = otherwise. 
This treatment reduces the melting temperature from in- 
finity to Tm- For this system we know from Eq. (^ that 
fj = 1/2, and we think it should be the same for the 
asymmetric model. 

Can the present approach be extended to consider the 
possible random variations in the on-site potential V{r) 
(this is caused by the sequence heterogeneity in the case 
of DNA jl^)? This is certainly a challenging problem 
and beyond the scope of this paper. But we think that 
inclusion of such an effect will not alter the qualitative 
behavior of the phase diagram, since in the renormaliza- 
tion sense, near the transition point the details of the 
interactions will be smoothed out |2^ . 

It is of interest to ask whether the inclusion of base- 
stacking effect (by making the parameter k in Eq. ([^) po- 
sition dependent as done in Ref. ) will change the ther- 
mal melting from second order to first order. It seems still 
be an issue of debate. An recent work done by Peyrard 
and coworkers pl[ | answered it confirmatively, while the 
numerical work of Cule and Hwa jl^ suggested that the 
transition behavior is still of second order. We noticed 
that in Ref. an force field is first included and at 
the final stage a limiting procedure is performed to make 
the field equal to zero. Our present work demonstrates 
that the property of the double-stranded system depends 
considerably on the external field, therefore it might be 
helpful for one to carefully evaluate whether the above 
mentioned limiting procedure in Ref. |2l[ causes a signif- 
icant effect. 

The author benefits from discussions with Xin Zhou 
and Yong Zhou. He is grateful to Z.-C. Ou-Yang for en- 
couragement and to W.-M. Zheng for bringing Ref. p9| ] 
to his attention. 
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FIG. 1. Phase diagram for a double-stranded biopoly- 
mer with a symmetric potentiaL The broken line 
shows a force-extension curve at temperature To. The 
phase-coexistence region collapses to a point only when 
T oo, indicating there is no thermal melting. (Inset) The 
relation between the average distance of the ends of the two 
strands and the force at temperature Tq. We set k = ksTo/b'^ 
and 7 = kBTodo- 




Monomer Extension a ((7/Ka) ) 

FIG. 2. Phase diagram for a double-stranded biopolymer 
with a asymmetric potential. This system shows second or- 
der thermal melting behavior at Tm ~ ^/nya/ks- Inset shows 
how the threshold force /c for melting changes with temper- 
ature T. 
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